Nonlocal bottleneck effect in two-dimensional turbulence 
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The bottleneck pileup in the energy spectrum is investigated for several two-dimensional (2D) 
turbulence systems by numerical simulation using high-order diffusion terms to amplify the effect, 
which is weak for normal diffusion. For 2D magnetohydrodynamic (MHD) turbulence, 2D electron 
MHD (EMHD) turbulence and 2D thermal convection, which all exhibit direct energy cascades, a 
nonlocal behavior is found resulting in a logarithmic enhancement of the spectrum. 

PACS: 47.27; 47.27Eq; 47.27Gs 



OO 

On 



3 



> ■ 

<n : 

o ■ 
■ 

o : 

00 . 
On ■ 

"a : 

■ 

i 1 

o 



x 



The local enhancement of the energy spectrum in front 
of the dissipation range, which is now generally called 
bottleneck effect, is a well-established phenomenon. It 
has been observed in numerous experiments [1],[2] and 
numerical simulations [3], [4], [5], and has been discussed 
theoretically [6] pointing out the physical mechanism. 
Even a quantitative formula was derived assuming a 
Batchelor fit for the second order structure function 
[7], [8]. It is interesting to note that the magnitude of 
the bottleneck effect, however, seems to depend on the 
character of the turbulent eddies. In simulations of super- 
sonic turbulence [9] the irrotational compressible part of 
the velocity field exhibits a considerably weaker spectral 
enhancement than the solenoidal part. For higher-order 
dissipation terms zA7 2 — > — v n {— V 2 )™, as often used in 
turbulence simulations to maximize the inertial range, 
the amplitude of the bottleneck effect increases, such 
that for n 3> 1 it seems to affect also the low-A: inertial 
range behavior [10], though these results are probably 
not asymptotic. 

Contrary to the attention the bottleneck effect at- 
tracted in three-dimensional (3D) turbulence, it has to 
our knowledge not yet been discussed in 2D turbulent 
systems. It is true that for the enstrophy cascade in 
2D Euler turbulence no such effect exists, which is at- 
tributed to the negative sign of the eddy viscosity [11], 
making the energy spectrum slightly steeper than the 
corresponding Kolmogorov law, Ek ~ fc~ 3 (ln fc/fco) -1 / 3 . 
But for 2D systems dominated by a direct energy cas- 
cade there is no a priori argument, why the same mech- 
anism leading to the bottleneck effect in 3D should not 
also be active in 2D. The effect seems, however, to be 
much weaker, since numerical simulations of such turbu- 
lent systems, in particular in 2D magnetohydrodynamics 
(MHD) [12], [13], [14] and 2D electron magnetohydrody- 
namics (EMHD) [15], [16] have found energy spectra ex- 
hibiting almost perfect power laws down to the dissipa- 
tive fall-off with no visible bottleneck pileup. 

One could argue that the difference is only due to the 
geometry of the triad interactions, now restricted to one 
plane, and that applying the analysis of Ref. [8] to a 2D 
system might lead to a less pronounced effect than in 3D. 



The only change in the algebra is to replace the integral 
in the expression for Ek oc kr sin krD(r)dr, see [8] , by 
/ °° krJo(kr)D(r)dr, where D(r) = ([u r (x + r) — v r (x)] 2 } 
is the second order longitudinal structure function. As- 
suming again a Batchelor fit for D(r), a straightforward 
evaluation shows that the bump on the spectrum would 
be of similar magnitude and width as in the 3D case, 
contrary to the much weaker effect revealed in 2D simu- 
lations, which invalidates this assumption. It is therefore 
necessary to investigate the character of the transition 
from the inertial to the dissipation range in 2D turbu- 
lence more in detail. 

In this letter we present results of a series of simula- 
tion runs for the turbulence systems mentioned above, 
2D MHD, 2D EMHD, and also 2D thermal convection 
[18], using high n in order to amplify the inherently weak 
2D bottleneck effect and going to higher spatial resolu- 
tion than done previously. All three systems are two- 
field models which, though formally of similar structure, 
exhibit rather different turbulence properties. Here the 
main interest is, however, not in the physics described 
by these models, for which we refer to the original pa- 
pers. We first consider EMHD turbulence, which is most 
closely related to (3D) Navier-Stokes turbulence. The 2D 
EMHD equations are [15] 



(d t + v e • V)(tf - d 2 e j) = - Vn (-V 2 r^, 



(d t + v e • V)(<£ - d 2 e cu) + B • Vj = -Vn(-V 2 



(1) 

(2) 



where the flux function ip describes the magnetic field 
in the plane, B = e z x Vip, j — V 2 V>, and the stream 
function (f> describes the electron flow in the plane, v e = 



x V</>, UJ 



v e is proportional to the current 



density in the plane, such that qb gives the out-of-plane 
field fluctuation, cj> = 8B Z . The equations are written in 
nondimcnsional form and d e — c/uj pe L is the normalized 
collisionless electron skin depth, for details see [15]. The 
equations are solved on a periodic box of linear size 27r 
using a standard pseudo-spectral method with dealias- 
ing according to the 2/3 rule. The dissipation terms 
are integrated exactly. As in [15] we consider turbu- 
lence decaying from a random initial state. It has been 
shown in [15] that for large wavenumbers kd e > 1, 2D 
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EMHD turbulence exhibits a Kolmogorov energy spec- 
trum E k ~ e 2 / 3 fc- 5 / 3 . 



is independent of the extent of the inertial range. Increas- 
ing the Reynolds number should only shift the bump to 
larger k, but not increase its amplitude, see e.g. [8]. We 




FIG. 1. Compensated energy spectra k 5/3 E k of 2D EMHD 
turbulence for diffusion operator order n = 3, 8, 20. Note the 
linear vertical scale. 

Figure |l| gives the energy spectrum E k = (E k e~ 2 / 3 ) t 
averaged over about one energy decay time. (Since at high 
n the dissipation length is essentially independent of e, 
the average can be performed at constant k [17].) Shown 
are three cases with d e = 0.3, iV 2 = 2048 2 , A: max = 682, 
chosen as in [15], (a) n — 3, 773 = 6 x 10~ n , (b) n — 
8, r; 8 = 1(T 38 , (c) n = 20, 7720 = 10" 96 . The modal en- 
ergy is defined by E k = E, ng i c (k 2 \^k\ 2 + \^\ 2 ){l+d 2 e k 2 ). 
While no bottleneck effect is visible for n — 3, in agree- 
ment with the spectrum shown in [15], there is a clear 
spectral enhancement for n — 8 of 20% and for n = 20 of 
about a factor of 2. For comparison with the correspond- 
ing 3D behavior several 3D EMHD simulation runs have 
been performed for similar parameter values, though at 
lower Reynolds numbers. 3D EMHD follows the equation 

d t (B - d 2 V 2 B) - V x [v e x (B - d 2 V 2 B)] 



-77 n (-V 2 ) n B, 



(3) 



where v e = —V x B. In Fig. || we plot the compensated 
energy spectrum from three simulation runs of decaying 
3D turbulence with N 3 = 256 3 , fc max = 85, for d e = 1 
and n — 3, 8, 20, which show bottleneck enhancement 
factors of 2.5,4,10, respectively. Hence the bottleneck 
effect is indeed quantitatively much weaker in 2D than 
in 3D. (Note that for kd e > 1 3D EMHD formally re- 
duces to the Navier-Stokes equation in the vorticity form 
V 2 B = V x v e — * V x v. In fact the n = 8 spectrum in 
Fig. H is practically identical with that observed for the 
corresponding Navier-Stokes case [17].) 

The general understanding of the bottleneck seems to 
be that the spectral enhancement, while depending on n, 
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FIG. 2. Compensated energy spectra E k for 3D EMHD 
turbulence simulations with n = 3, 8, 20. Normalization is 
such that the horizontal parts of the spectra coincide. 



will now show that such behavior is in general not true 
for 2D turbulence. Figure || gives the energy spectrum 
for three 2D EMHD turbulence runs, where we choose 
n = 8 to amplify the effect, ry 8 = 10~ 32 , 10~ 36 , 10~ 40 us- 
ing N 2 = 1024 2 ,2048 2 ,4096 2 , respectively. While there 
is no visible bottleneck effect for r]s — 10~ 32 , it becomes 
more and more pronounced with decreasing dissipation 
coefficient leading to a flattening of the spectrum E k 
in an increasingly larger fraction of the inertial range, 
k > k ~ 60. 




FIG. 3. Compensated energy spectra for three 2D EMHD 
turbulence simulations with n = 8, r/ 8 = 10~ 32 , 10~ 36 , 10~ 40 . 

This behavior is not limited to 2D EMHD turbulence, 
but is found to occur in a similar and even clearer form 
in 2D MHD turbulence simulations. Here the dynamical 
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equations are 

d t ip + v • V^> 
d t LO + v ■ Vw — B • Vj 



-77„(-V 2 )>, (4) 
-^(-V 2 )^, (5) 



where v = e z x is the plasma flow (note that 
in spite of the formal similarity EMHD does not con- 
verge to MHD for d e < 1. EMHD is limited to 
d e > ^Jm e /mi, while MHD is only valid at larger, 
macroscopic scales, see [15]). Previous numerical stud- 
ies of decaying 2D MHD turbulence have revealed the 
spectral law E k ~ (w^e) 1 / 2 ^^ 3 / 2 (see e.g. [14]), where 
va = B/y/Anp is the Alfven speed and the modal en- 
ergy is E k = £ anglc k 2 {\^ k | 2 + \<f> k \ 2 ). The fc" 3 / 2 spec- 
trum results from the Alfven effect [19], the coupling of 
small-scale velocity and magnetic field fluctuations by 
the magnetic field of the large-scale eddies. For normal 
diffusion no bottleneck is discernable in the energy spec- 
tra [12], [14]. To investigate this point more closely, we 
choose again a high-order dissipation operator in order 
to amplify the bottleneck effect, which may be hidden in 
the noise level for n = 1. Three simulation runs have 
been performed for decaying MHD turbulence using the 
same numerical scheme as in the EMHD simulations de- 
scribed above and similar initial conditions (called B-type 
in [12]). Figure | gives the time-averaged MHD energy 
spectra E k = (E k (vA£)~ 1 ^ 2 )t plotted in compensated 
form for rj 8 = i/$ = 10" 36 , 10~ 40 , 10~ 45 with resolutions 
N 2 = 1024 2 ,2048 2 ,4096 2 , respectively. Comparing with 
the corresponding 2D EMHD cases given in Fig. 0, Fig. 
^ shows a similar qualitative trend, a nonlocal influence 
of the dissipation range on the inertial range. The differ- 
ence is probably due to the choice d e — 0.3 in the EMHD 
runs in Fig. ^. Since a pure scaling behavior exists only 
for kd e > 1, the effective scaling range is shorter by a 
factor 3-4 for the same resolution, such that amplitude 
of the bottleneck pileup in the highest resolution EMHD 
case 4096 2 in Fig. || corresponds to the lowest-resolution 
MHD case 1024 2 in Fig. g 

There is a nearly linear increase of the compensated 
spectrum in the log-linear plot, which suggests that for 
sufficiently large k > k ~ 20 the inertial range is modified 
by a logarithmic factor 

£ fc ~£T 3 / 2 ln(fc/ K ), (6) 

where the magnitude of the effect is expected to depend 
on n, such that for n = 1 it becomes invisibly small at the 
achievable spatial resolution. The possibility of a loga- 
rithmic factor has been discussed in [6] for 3D turbulence, 
though only as a subdominant effect in the spectral cor- 
rection term. The wavenumber k is connected with some 
structure of the macro-state of the system. The fact that 
k is lower in the MHD runs than observed in the EMHD 
runs is due to the choice d e = 0.3 in the latter. 

These results show that the bottleneck effect in the 2D 
turbulence systems considered, though very weak for nor- 
mal diffusion n = 1, exhibits a nonlocal behavior when 




FIG. 4. Compensated 2D MHD energy spectra k 3/2 E k for 



r ?8 = ^ 8 = 10 - 36 ,10- 40 ,10' 45 . 



enhanced by choosing high n. The mechanism for this 
property must be connected with a stronger direct in- 
teraction of small- and large-scale modes in 2D than in 
3D, which is also reflected in configuration space by the 
large-scale intermittency typical for 2D turbulence. The 
observed behavior is probably not due to the Alfven ef- 
fect, since the latter is not present in EMHD [15]. 

We would like to discuss briefly also a third type of 
turbulence, 2D thermal convection in the Boussinesq ap- 
proximation described by the equations 



StT + vVTl d v <t> = -x„(-V 2 ) n T, 



dtw + v • Vw 



d y T = 



(7) 
(8) 



again written in nondimensional form, see [18], whre T 
is the temperature fluctuation and v — e z x V</> the fluid 
velocity. Contrary to MHD or EMHD, where stationary 
turbulence can only be achieved by an external stirring 
force, this system is linearly unstable over a broad k- 
range with growth rate 7 oc k y /k, which generates a sta- 
tionary level of turbulence. This is caused by the frozen- 
in mean temperature gradient Tg, v x Tq = d v 4> in the 
nondimensional form (|t]). Hence there is no ideal energy 
invariant, instead one has 



v 2 )d 2 x = 2 



r T d 2 x 



(9) 



where e is the energy dissipation rate. Equations (0),(||) 
have recently been studied numerically on a periodic 
box using a similar scheme as in the EMHD simula- 
tions. In spite of the anisotropic linear drive the spec- 
tra 2?J = \T k \ 2 and E\ — |i; fe | 2 are highly isotropic, 
which demonstrates the strong influence of the nonlinear 
terms. The energy dissipation rate e — + + ejf + 
in (||) consists of the dissipation on the velocity and the 
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temperature fluctuations both at large (L) and small (s) 
wavenumbers. The kinetic energy v 2 has an inverse cas- 
cade and is primarily dissipated at small wavenumbers, 
where the modes are artificially damped to prevent con- 
densation and suppression of turbulence [18], < e]f , 
while the thermal fluctuation energy T 2 has a direct cas- 
cade, > ej . 

The results of [18] seem to indicate spectral laws 
E\ ~ fc~ 12 and E\ ~ k~ 2 - 3 , which differ somewhat 
from the expected Bolgiano scaling [20] k~ 1A and fc~ 2,2 , 
respectively. No convincing physical argument for these 
deviations could be given in [18] except the fact that the 
temperature fluctuations are found to be highly intermit- 
tent, which limits the relevance of the (nonintermittent) 
Bolgiano scalings SiT ~ I - 2 , Siv ~ I - 6 . Here we suggest 
an alternative interpretation of the simulation results in 
[18], namely that the deviations are caused by nonlocal 
effects, which should make the spectrum of the kinetic 
energy with an inverse cascade slightly steeper (cf. the 
2D Euler case), while the spectrum of the temperature 
fluctuations, which have a direct cascade, should be flat- 
ter than the Bolgiano spectrum. For direct comparison 
with the EMHD and MHD simulations presented above, 
we have performed a similar run for thermal convection 
with n = 8, xs = v % = 10~ 42 , N 2 = 2048 2 , from which 
Fig. ^ shows the temperature fluctuation spectrum com- 
pensated with the Bolgiano law. The behavior is indeed 
very similar to the MHD spectrum in Fig. [|. 

In conclusion we have investigated the spectral bot- 
tleneck pileup in 2D turbulence. While it is known 
that there is no such effect in the enstrophy cascade 
for 2D Euler turbulence, we have shown the existence 
of the effect in 2D systems exhibiting a direct energy 
cascade, in particular in MHD, EMHD and thermal con- 
vection. The amplitude of the pileup is found to be sig- 
nificantly smaller than in corresponding 3D cases even for 
high-order hyperdiffusion, which enforces the tendency of 
spectral pileup by making the transition between inertial 
range and dissipative range more abrupt. But contrary 
to the local behavior in 3D the bottleneck effect in 2D 
turbulence has a nonlocal character, such that the ma- 
jor part of the inertial range is affected for sufficiently 
small dissipation coefficients. The simulation results sug- 
gest a logarithmic modification of the spectrum. This 
behavior reflects a strong direct interaction of disparate 
scales, which also manifests itself in the large-scale in- 
termittency typically seen in configuration space plots of 
2D turbulence. Also the deviations in the spectral laws 
from Bolgiano scaling in 2D thermal convection reported 
previously can be attributed to this nonlocal bottleneck 
effect. As a consequence the popular use of high-order 
diffusion operators in 2D turbulence simulations becomes 
rather doubtful. 

The authors would like to thank Andreas Zeiler for 
providing the basic version of the 3D EMHD code. A.C. 
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FIG. 5. 2D thermal convection. Compensated temperature 
fluctuation spectrum k ' for n — 8. 
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